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C^j Results are presented on event-by-event fluctuations of transverse momenta pr in 

central Pb+Pb interactions at 20 A, 30A, 40A, 80A, and 158A GeV. The analysis was 
performed for charged particles at forward center-of-mass rapidity (1.1 < y* < 2.6). 
Three fluctuation measures were studied: the distribution of average transverse 
momentum - M(pr) - in the event, the & PT fluctuation measure and two-particle 
transverse momentum correlations. Fluctuations of px are small and show no signifi- 
cant energy dependence in the energy range of the CERN Super Proton Synchrotron. 
Results are compared with QCD-inspired predictions for the critical point, and with 
the UrQMD model. Transverse momentum fluctuations, similar to multiplicity fluc- 
tuations, do not show the increase expected for freeze-out near the critical point of 
QCD. 
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1 Introduction and Motivation 



For more than 30 years, experiments studying relativistic nucleus-nucleus (A + A) colli- 
sions have been carried out in laboratories in Europe and the United States. The main 
motivation has been to test the hypothesis that strongly interacting matter at energy 
densities exceeding about 1 GeV/fm 3 exists in the form of deconfmed quarks and gluons, 
eventually forming the quark-gluon plasma (QGP) [1]. Recent results from the CERN 
Super Proton Synchrotron (SPS) and BNL Relativistic Heavy Ion Collider (RHIC) seem 
to confirm this conjecture. The data suggest that the threshold for the onset of decon- 
finement is located at the low SPS energies [2, 3]. 

The phase diagram of strongly interacting matter can be presented in terms of tem- 
perature T and baryochemical potential hb- QCD-inspired calculations suggest that the 
phase boundary between hadrons and QGP is of first order at large values of /j,b, ending 
in a critical point of second order and then turning into a continuous rapid transition at 
low (j,b [4]. The location of the critical point may be close to the (T,/ib) values found for 
the freeze-out of the hadron system produced in collisions of heavy nuclei at SPS energies. 

The study of fluctuations is an important tool for localizing the phase boundary and 
the critical point. In particular, significant transverse momentum and multiplicity fluc- 
tuations are expected to appear for systems that hadronize from QGP and freeze-out 
near the critical point of QCD [5]. The location of the freeze-out point in the phase di- 
agram can be moved by varying the collision energy and the size of the collision system. 
A nonmonotonic evolution of fluctuations with these collision parameters can serve as a 
signature for the phase transition and the critical point. These considerations motivated 
an extensive program of fluctuation studies at the SPS and RHIC accelerators. 

So far for central Pb+Pb collisions, the analysis of multiplicity fluctuations has found 
only small effects without significant structure in the energy dependence in the whole SPS 
energy domain [6]. The measured net charge fluctuations can be explained by the effects 
of global charge conservation [7]. The energy dependence of event-by-event fluctuations 
of the K/ir ratio, on the other hand, shows an interesting increase toward lower energies 
[8] that can be attributed to the onset of deconfinement [9] rather than to the critical 
point. 

This paper presents results of the NA49 experiment from a study of transverse mo- 
mentum pt fluctuations in central Pb+Pb collisions at 20A, 30A, 40A, 80A, and 158A 
GeV. It extends a previous study [10] that investigated the system size dependence of pr 
fluctuations at the top SPS energy. 

Fluctuations in nucleus-nucleus collisions are susceptible to two trivial sources: the 
finite and fluctuating number of produced particles and event-by-event fluctuations of the 
collision geometry (see the discussion of this point within HSD and UrQMD transport 
models in Ref. [11]). Suitable statistical tools have to be chosen to extract the fluctuations 
of interest. As in the previous NA49 study [10], mainly the global $ PT measure will be 
used. Alternative measures used in the literature, e.g. o VT ^ yn [12], Aa PT [13], F PT [14], 
and E PT (%) [15, 16] are related [17]. 

To obtain further information on the possible source of fluctuations, two additional 
methods will be employed. The distribution of event-wise average transverse momentum 
M{pt) will be compared against the corresponding histogram for mixed events, which 
represents purely statistical fluctuations. Moreover, two-particle transverse momentum 
correlations, as proposed in Ref. [13], will be analyzed. 

This paper is organized as follows. In Sec. 2 the statistical tools used in this analysis 
are introduced and briefly discussed. The NA49 setup is presented in Sec. 3. Experimen- 
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tal effects such as detector acceptance and two-track resolution are discussed in Sec. 4. 
The NA49 results on the energy dependence of transverse momentum fluctuations are 
presented and discussed in Sec. 5. A summary closes the paper. 



2 Measures of fluctuations 

Various methods can be used to measure event-by-event px fluctuations. A natural ob- 
servable is the distribution of the average transverse momentum of the events defined 
as 

1 N 

M(pr) = TT^PTi, (1) 
ly i=l 

where N is the multiplicity of accepted particles in a given event and pn is the transverse 
momentum of the i-th particle. The distribution of M{p T ) will be compared with the 
corresponding histogram obtained for artificially produced "mixed events". In mixed 
events, all particles are by construction uncorrelated (each particle in a given mixed event 
is taken from a different real event with the same multiplicity) but follow the experimental 
inclusive three-momentum spectra as well as the distribution of multiplicity. 

The second observable used in this work is the $ PT measure, proposed in Ref. [18] 
and used also in our previous analysis [10]. Following the authors of Ref. [18], one can 
define the single-particle variable z PT = pr ~Pt with the bar denoting averaging over the 
single-particle inclusive distribution. As = 0. Further, one introduces the event 

variable Z PT , which is a multiparticle analog of z PT , defined as 

N 

Z pt = Z!(Pri -Pt), (2) 
1=1 

where the summation runs over particles in a given event. Note, that (Z PT ) = 0, where 
(...) represents averaging over events. Finally, the $ PT measure is defined as 



(Z PT ) /== 



\ (AO 



(3) 



The $ PT measure has two important properties. First, <3> PT vanishes when the system 
consists of particles that are emitted independently (no interparticle correlations), and 
the single particle momentum spectrum is independent of multiplicity. Second, if an 
A + A collision can be treated as an incoherent superposition of independent N + N 
interactions (superposition model), then <3> PT has a constant value, the same for A + 
A and N + N interactions. This implies that, in particular, $ PT does not depend on 
the impact parameter (centrality), if the A + A collision is a simple superposition of 
N + N interactions. Furthermore, $ PT is independent of changes of the size of a compact 
acceptance domain provided the correlation scale (range), lc, is much smaller than the size 
of the acceptance region, I a- In the limit of large correlation length and small acceptance 
the magnitude of $ PT is proportional to the multiplicity of accepted particles. Thus for a 
small acceptance $ PT approaches zero. The approximately linear dependence of $ PT on 
the fraction of accepted particles in the limit I a lc suggested the introduction of the 
fluctuation measures E PT (%) [16] and <J PT ^ yn [12], which for large particle multiplicities 



are proportional to JQ PT /(N). These properties of the fluctuation measures should be 
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taken into account when results are compared and discussed. Finally, we note that the 
magnitude of $ PT decreases with the fraction of randomly lost particles, e.g., due to 
incomplete reconstruction efficiency 

In spite of the above mentioned advantages, there is an important disadvantage of 
using <£> PT in the fluctuation analysis. While <3> PT is sensitive to the presence of particle 
correlations in a system, it does not provide information on the nature of the correlation. 
Several effects can contribute to <3> PT . Therefore, to achieve a better understanding of the 
fluctuation structure, it is useful to also employ a more differential method [13]. 

The correlations can be studied by plotting the cumulative pr variables for particle 
pairs. Namely, instead of pr, one can introduce the variable x, defined for a particle i as 



where p(pr) is the inclusive pt distribution, normalized to unity, which is obtained from 
all particles used in the analysis. By construction, the x variable varies between and 
1 with a flat probability distribution. Two-particle correlation plots, as presented in this 
paper, are obtained by plotting (xi,x 2 ) points for all possible particle pairs within the 
same event. The number of pairs in each (x 1 ,x 2 ) bin is divided by the mean number of 
pairs in a bin (averaged over all (xi,x 2 ) bins). This two-dimensional plot is uniformly 
populated when no interparticle correlations are present in the system. Nonuniformity 
signals the presence of interparticle correlations. For example, when identical particles are 
studied, Bose statistics lead to a ridge along the diagonal of the (x±, x 2 ) plot, which starts 
at (0,0) and ends at (1,1), whereas event-by-event temperature fluctuations produce a 
saddle-shaped structure [13]. 

3 Experimental Setup 

NA49 is one of the fixed target experiments at the CERN SPS. The detector (see Fig. 1 
and Ref. [20]) is a large acceptance hadron spectrometer used for the study of the hadronic 
final states produced in p + p, p+nucleus, and nucleus+nucleus collisions. In particular, 
centrality selected Pb+Pb interactions were recorded at 20A, 30A, AOA, 80A, and 158A 
GeV projectile energies. 

The main components of the detector are four large- volume time projection chambers 
(TPCs), which are capable of detecting about 80% of approximately 1500 charged particles 
created in a central Pb+Pb collision at 158A GeV (the acceptance losses are concentrated 
mainly in the backward rapidity region and at azimuthal angles along the magnetic field 
direction). The vertex TPCs (VTPC-1 and VTPC-2), are located in the magnetic field of 
two superconducting dipole magnets (1.5 and 1.1 T, respectively, at 158A GeV; for lower 
energies, the magnetic field is scaled down in proportion to the beam energy). Two other 
TPCs (MTPC-L and MTPC-R) are positioned downstream of the magnets symmetrically 
to the beam line. The results presented in this paper are analyzed with a global tracking 
scheme [21], which combines track segments that belong to the same physical particle 
but were detected in different TPCs. The NA49 TPCs allow precise measurements of 
particle momenta p with a resolution of o~(p)/p 2 = (0.3 — 7) • 10~ 4 (GeV/c) -1 . A precise 
measurement of specific energy loss (dE/dx) in the region of relativistic rise is possible in 
the TPCs, however, dE/dx information is not used in this analysis. 

The centrality of nuclear collisions is selected by using the energy of the projectile 
spectator nucleons measured in the downstream calorimeter (VCAL). The geometrical 
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Figure 1: Experimental setup of the NA49 experiment [20] with different beam definitions and 
target arrangements. 

acceptance of the VCAL is adjusted to cover the projectile spectator region by setting 
the collimator (COLL). 

The target is a thin Pb foil (224 mg/cm 2 ) positioned about 80 cm upstream from 
VTPC-1. Pb beam particles are identified by means of their charge as seen by a counter 
(S2') situated in front of the target. The beam position detectors (BPD-1/2/3 in Fig. 1), 
which are proportional chambers placed along the beam line, provide a precise measure- 
ment of the transverse positions of the incoming beam particles. For Pb beams, interac- 
tions in the target are selected by an anti-coincidence of the incoming beam particle with 
a counter (S3) placed directly behind the target. 

Details of the NA49 detector setup and performance of the tracking software are 
described in [20]. 

4 Data selection and analysis 
4.1 Data sets 

The data used for the analysis consist of samples of Pb+Pb collisions at 20A, 30A, 40A, 
80A, and 158A GeV energy for which the 7.2% most central reactions were selected. 
Table 1 describes the data sets used in this analysis. 



Collision energy (GeV) 


Year 


No. events 


y/s NN (GeV) 


Vera. 


20 A 


2002 


229 000 


6.27 


1.88 


30 A 


2002 


297 000 


7.62 


2.08 


40 A 


1999 


165 000 


8.73 


2.22 


80 A 


2000 


228 000 


12.3 


2.56 


158 A 


1996 


166 000 


17.3 


2.91 



Table 1: Data sets used in this analysis, for 7.2% central Pb+Pb collisions: collision energy, year 
of data taking, number of events, center-of-mass energy y / S7v r /v for N+ N pair, and center-of-mass 
rapidity y c . m . in the laboratory frame. 
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4.2 Event and particle selection 



Event selection criteria were aimed at reducing possible contamination with nontarget 
collisions. The primary vertex was reconstructed by fitting the intersection point of the 
measured particle trajectories. Only events with a proper quality and position of the 
reconstructed vertex are accepted in this analysis. The vertex coordinate z along the 
beam has to satisfy \z — z \ < Az, where zo is the nominal vertex position and Az is a cut 
parameter. The values of z and Az are -581.05 and 0.25 cm for 20A GeV, -581.3 and 0.3 
cm for 30A GeV, -581.1 and 0.3 cm for AO A GeV, -581.2 and 0.3 cm for 80A GeV, and 
-578.9 and 0.3 cm for 158A GeV. The maximal allowed deviation of the x and y positions 
of the fitted vertex (Ax and Ay) from the measured beam position varies from 0.1 cm (for 
158A GeV) to 0.25 cm (for low energies). The cut on the ratio of the number of tracks 
used to fit the primary vertex divided by the total number of tracks registered in TPCs 
(ntf /nto) was required to be higher than 0.25 for all energies. 

To reduce the possible contamination by nonvertex tracks such as particles from sec- 
ondary interactions or from weak decays, several track cuts are applied. The accepted 
particles are required to have measured points in at least one of the vertex TPCs. A cut 
on the so-called track impact parameter, the distance between the reconstructed main 
vertex and the back extrapolated track in the target plane \b x \ < 2 cm and \b y \ < 1 cm, is 
applied to reduce the contribution of nonvertex particles. Moreover, particles are accepted 
only when the potential number of points (nmp), calculated on the basis of the geometry 
of the track, in the detector exceeds 30. The ratio of the number of points on a track to 
the potential number of points (np/nmp) is required to be higher than 0.5 to avoid using 
fragments of tracks. This set of cuts significantly reduces the possible contamination of 
particles from weak decays, secondary interactions, and other sources. 

In this analysis only tracks with 0.005 < pt < 1.5 GeV/c are used. For all five energies, 
the forward-rapidity region is selected as 1.1 < < 2.6, where y* is the particle rapidity 
calculated in the center-of-mass reference system. As the track-by-track identification is 
not always possible in the experiment, the rapidities are calculated assuming the pion 
mass for all particles. 

The NA49 detector provides a large, though incomplete, acceptance in the forward 
hemisphere. Figure 2 presents examples of ((f), Pt) acceptance 1 for 2.0 < y* < 2.2. 
The regions of complete azimuthal acceptance are different for various energies. The 
quantitative comparison of $ PT values among the five energies requires selection of the 
regions of common acceptance. 

The solid lines in Fig. 2 represent a parametrization of the common acceptance limits 
for all five energies by the formula - 

M0) = j 2 - B, (5) 

where the parameters A and B depend on the rapidity range as given in Table 2. The 
same limits are used for both the negatively charged and positively charged particles with 
the appropriate redefinition of azimuthal angle 1 . Only particles within the acceptance 
limits are used in the analysis. 

1 All charged particles are plotted, with the azimuthal angle of negatively charged particles (assuming 
standard polarity of the magnetic field) reflected: namely for all negatively charged particles with <j> < 
degrees their azimuthal angle is changed as follows: <f> goes to 0+360 degrees, and finally for all negatively 
charged particles cj> goes to <fi — 180 degrees. In the case of the opposite polarity of the magnetic field 
(40A GeV data) one has to redefine the azimuthal angle of positively charged particles, instead. 
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Figure 2: (Color online) NA49 {4>,Vt) acceptance of all charged particles for 2.0 < y* < 2.2. 
Additional cut on y* (see the text) is not included. The solid lines represent the analytical 
parametrization of the common acceptance. 





1.0-1.2 


1.2-1.4 


1.4-1.6 


1.6-1.8 


1.8-2.0 


2.0-2.2 


2.2-2.4 


2.4-2.6 




600 


700 


1000 


2600 


3000 


2500 


1800 


1000 




0.2 


0.2 


0.2 


0.5 


0.4 


0.3 


0.3 


0.1 



Table 2: Parametrization of the NA49 <j) — pt acceptance common for all five energies for 
positively charged particles (standard configuration of magnetic field). For negatively charged 
particles one has to change the definition of the azimuthal angle (see text) and then use the 
same parametrization. 



Due to the described kinematic cuts (mostly on rapidity) and track selection criteria, 
together with the requirement to use the same limited (4>,Pt) acceptance, only about 5.0 
- 5.2 % of charged particles produced in central Pb+Pb interactions (at each of the five 
energies) enter into the subsequent analysis. 

Figure 3 presents (y*,pr) plots of all charged particles accepted in the analysis (addi- 
tional cut on y* - see below - not applied in the plots). Top, middle and bottom panels 
correspond to 20A, 30A, and 80A GeV data, respectively, whereas left, middle, and right 
panels are obtained by assuming pion, kaon and proton masses when evaluating rapidities. 
Black lines represent beam rapidities (yl eam ) in the center-of-mass reference system. At 
lower energies the NA49 TPC acceptance extends to the projectile rapidity domain and 
the selected particles may be contaminated by e.g. elastically scattered or diffractively 
produced protons. This domain has been excluded from the analysis by applying an ad- 
ditional cut on the rapidity y* of the particles, calculated with the proton mass. Namely, 
at each energy, the rapidity y* is required to be lower than yl eam — 0.5. 

The above cut was introduced because the preliminary analysis (without the additional 
y* cut) manifested an unexpected effect: the $ PT measure showed a different behavior for 
particles of different charges [22]. $ PT was independent of energy and consistent with zero 
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accept, part, at 20A GeV 




accept, part, at 80A GeV 




Figure 3: (Color online) (y*,pr) plots of all accepted (see Fig. 2) particles assuming pion (left), 
kaon (middle) and proton (right) mass. Additional cut on y* (see the text) is not included. Top, 
middle, and bottom panels correspond to 20^4, 30^4, and 80A GeV data, respectively. Black 
lines represent beam rapidities {yl eam ) in the center-of-mass reference system. 

for negatively charged particles, but it significantly increased for lower SPS energies for 
both all charged and positively charged particles. This effect, however, was found to be 
connected with protons only, and can be explained by event-by-event impact parameter 
fluctuations or, more precisely, by a correlation between the number of protons (nucleons) 
in the forward hemisphere and the number of protons (nucleons) that are closer to the 
production region [23] . For more central events the number of forward-rapidity protons is 
smaller, and consequently, the number of protons in the production region is higher. The 
situation is opposite for less central collisions. The existence of those different event classes 
results in the increased $ PT values for positively charged particles. One can eliminate 
this trivial source of correlations by either centrality restriction or rejection of the beam 
rapidity region [23]. In this analysis, the second method is employed by applying a cut 
on y* at each energy (see above). 
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4.3 Corrections and error estimates 

The statistical error on $ PT has been estimated as follows. The whole event sample 
was divided into 30 independent subsamples. The value of $ PT was evaluated for each 
subsample, and the dispersion - D - of the results was then calculated. The statistical 
error of $ PT is taken to be equal to D / ^/30. 

The event and track selection criteria reduce the possible systematic bias of the mea- 
sured <3> PT values. To estimate the remaining systematic uncertainty the values of cut 
parameters were varied within a reasonable range. For a given cut parameter, the "par- 
tial" systematic error was taken as half of the difference between the highest and the 
lowest $ PT value. Two event cuts and two track cuts were considered in the analysis. 
The final estimate of the total systematic error on $ PT was taken as the maximum of the 
changes resulting from this study. 

Event cuts are used to reject possible contamination of nontarget interactions, however 
there is always a small fraction of remaining nontarget events that can influence the $ PT 
values. This systematic bias can be estimated by investigating the dependence of $ PT on 
two vertex cut parameters - ntf /nto and Az. Figure 4 presents the dependence of $ PT on 
the ratio ntf /nto for 20A, 30 A, and 80A GeV data. The observed systematic error of $ PT 
with respect to changes of ntf /nto varies from 0.35 MeV/c for 20A GeV to 0.85 MeV/c 
for 30A GeV and 80A GeV. In Fig. 5 the dependence of $ PT on the allowed distance Az 
from the nominal position of the main vertex is shown for three energies. The observed 
variation of $ PT with this cut is quite small for all studied data sets. The highest "partial" 
systematic error (for 80A GeV data) was found to be 0.75 MeV/c. 
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Figure 4: Dependence of <& PT on one of the event cuts (ntf/nto ratio). Note: the values and 
their errors are correlated. The dashed line indicates the cut used in the analysis. 



The majority of tracks selected by the track selection criteria are main vertex tracks 
and the remaining fraction (?al0%) originates predominantly from weak decays and sec- 
ondary interactions with the material of the detector. The influence of this remaining 
fraction can be estimated by studying the dependence of $ PT on the track cut parameters 
\b x \ and \b y \ and the np/nmp ratio. Figures 6 and 7 present how the values of $ PT change 
with the impact parameter cut and with the cut on the np/nmp ratio, respectively. One 
finds that <3> PT is rather stable with respect to both cut parameters. A small increase of 
$ PT with increasing impact parameter cut may be due to the increasing contribution of 
nonvertex tracks for higher \b x \ and \b y \ values. This effect has been studied quantitatively, 
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• 20A GeV 
■ 30A GeV 

' 80A GeV 



j .... . 

T x ! 1 1 1 1 1 1 1 1 



0.1 0.2 0^3 0.4 0^5 0^6 0^7 

Az [cm] 

Figure 5: Dependence of <& PT on the allowed distance Az from the nominal position of the 
main vertex. Note: the values and their errors are correlated. The dashed line indicates the cut 
used in the analysis of 30A and 80A GeV data. The dotted line represents the cut used in the 
analysis of 20A GeV data. 

for central Pb+Pb collisions at 158 A GeV, in our previous paper [10]. The systematic 
error contribution from the impact parameter cut dependence varies from 0.15 MeV/c for 
80 A GeV to 0.55 MeV/c for 20 A GeV, whereas the contribution from the change of cut in 
the np/nmp ratio varies from 0.45 MeV/c for 40A GeV and 80A GeV up to 1.45 MeV/c 
for 158A GeV. 




• 20A GeV 
■ 30A GeV 

' 80A GeV 

I I T T 1 1 I T T 

ft " i i i i 



1 1.5 2 2.5 3 3.5 4 4.5 
cut on |bj [cm] 

Figure 6: Dependence of <3? PT on upper cut in the impact parameter \b x \. For each point, the 
cut on \b y \ is equal to half the cut on \b x \. Note: the values and their errors are correlated. The 
dashed line indicates the cut used in the analysis. 

Finally, the typical systematic error on $ PT , determined as the maximum resulting 
from the above analysis, has been estimated as not higher than 0.7 MeV/c for 20A GeV, 
0.9 MeV/c for 80A GeV, 1.1 MeV/c for 40A GeV, and 1.6 MeV/c and 1.7 MeV/c for 30A 
and 158A GeV data, respectively. More detailed values (also for various charge selections) 
are given in Table 4. 

The NA49 experiment registered 158A GeV Pb+Pb interactions at various beam in- 
tensities, magnetic field configurations ("normal" polarity of the magnetic field STD+ 




11 



> 

ID 



T3 
O 
U> 
i- 

.c 

o 



• 20AGeV 
■ 30AGeV 
' 80AGeV 



0.2 0.4 0.6 0.8 1 1.2 
cut on np/nmp 



Figure 7: Dependence of & PT on lower cut on np/nmp ratio. Note: the values and their errors 
are correlated. The dashed line indicates the cut used in the analysis. 



and opposite polarity of the magnetic field STD-), and with different centrality triggers. 
Figure 8 presents <3> PT values for these different data sets. An additional cut on the energy 
deposited in VCAL has been applied to select the 7.2% most central Pb+Pb interactions 
from each data sample. Although the numbers of events extracted from minimum-bias 
data are relatively low, one observes good agreement between the Q PT values obtained for 
different experimental conditions. 
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Figure 8: 3> PT values for different central and minimum-bias (MB) data sets. Data used for the 
further analysis are central STD+ normal intensity. 



The observed <£> PT values can be affected by the losses of tracks resulting from the 
reconstruction inefficiency of the detector and from the track selection cuts. The depen- 
dence of $ PT on the fraction of randomly rejected particles is shown in Fig. 9. Within 
the considered forward-rapidity region, at all five energies (only three example energies 
are shown), the tracking efficiency of the NA49 detector is higher than 95% and therefore 
Fig. 9 implies that the bias due to tracking inefficiency and track selection cuts is lower 
than 0.5 - 1.0 MeV/c. 

It has already been shown in our previous paper [10] that the limited two-track reso- 
lution influences the measured <3> PT values. To quantitatively estimate this contribution, 
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Figure 9: Dependence of <& PT on the fraction of randomly rejected particles. 

five samples (30 000 events in each) of mixed events were prepared (for 20A, 30A, 40A, 
80 A, and 158 A GeV data). Multiplicities of mixed events were chosen to be the same 
as those of real events, but each particle in a mixed event was taken at random from 
a different real event. Such a mixing procedure ensures that the inclusive spectra (e.g., 
transverse momentum or rapidity) are the same for data and for mixed events, but all 
possible correlations vanish. Indeed, it was verified for all data sets and all particle charge 
selections that the $ PT value calculated for the sample of mixed events is consistent with 
zero. In a second step, the prepared mixed events were processed by the NA49 simulation 
software. The resulting simulated raw data were then reconstructed and the $ PT measure 
calculated. The obtained $ PT values are negative, as expected for the anticorrelation 
introduced by the losses due to the limited two-track resolution. The additive two-track 
resolution correction is calculated as the difference (A$ PT ) between the values of $ PT 
after detector simulation and reconstruction and before this procedure (mixed events). 
The values of A<3> PT were calculated for all five energies and for various particle charge 
selections, separately, and are listed in Table 3. 



Particles 


20A 


30A 


A0A 


80A 


158A 


All 


-0.6 ± 0.4 


-0.1 ± 0.3 


-0.9 ± 0.4 


-1.7 ± 0.4 


-3.4 ± 0.6 


Negative 


-0.1 ± 0.2 


0.4 ± 0.2 


-0.7 ± 0.3 


-1.4 ± 0.2 


-1.7 ± 0.5 


Positive 


-0.2 ± 0.4 


-0.5 ± 0.4 


-0.6 ± 0.4 


-1.8 ± 0.5 


-4.2 ± 0.7 



Table 3: Values of two-track resolution corrections A<E> PT (in MeV/c) for all accepted charged 
particles and separately for negatively and positively charged particles, for each of the five 
energies (in GeV). 

The magnitude of the two-track resolution corrections is higher for higher energies 
(from 2 to 4 MeV/c for 158A GeV) where the multiplicities of produced particles are higher 
and consequently the densities of tracks are relatively high. The two-track resolution 
corrections become much smaller for low-multiplicity interactions, namely for 20A and 
30 A GeV data (A$ PT values close to zero). The absolute values of A$ PT are typically 
higher for positively than for negatively charged particles, which is mainly due to higher 
track density for positively charged particles caused by the larger number of protons 
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relative to antiprotons. The A$ PT values are typically negative thus indicating that <3? PT 
measured with an ideal detector would be higher. The value of $ PT corrected for the 
limited two-track resolution effect equals the "raw" $ PT minus the corresponding A$ PT . 
The error of the corrected <E> PT is calculated by adding in squares the statistical error of 
the raw <3> PT value and the statistical error of the correction. As A$ PT is very small for 
20A and 30A GeV data, no corrections for the two-track resolution effect were applied at 
these energies. Instead, the upper limits of the systematic errors were increased by 0.3 
MeV/c (see Table 4). 



5 Results and Discussion 
5.1 Results 

The results shown in this section refer to accepted particles, i.e., particles that are accepted 
by the detector and pass all kinematic cuts and track selection criteria as discussed in 
Sec. 4.2. Results are not corrected for limited kinematic acceptance, and this acceptance 
has to be taken into account when the data are compared with model predictions. The 
measured $ PT values are corrected for limited two-track resolution of the NA49 detector 
(see below and Sec. 4.3). A possible bias due to particle losses and contamination in the 
accepted kinematic region is included in the systematic error. 
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Figure 10: (Color online) Distributions of mean transverse momentum for data (points) and 
mixed events (histograms). Data points are not corrected for acceptance and limited two-track 
resolution. Events with accepted particle multiplicity equal to zero are not used. 



Figure 10 shows the distributions of mean (per event) transverse momentum M(px) 
for central 20A, 30A, 40A, 80A and 158A GeV Pb+Pb interactions. Points represent 
results from real events, and histograms are calculated from mixed events, in which there 
are no interparticle correlations by construction. Data are not corrected for experimental 
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effects such as two-track resolution or acceptance. Events with zero multiplicity (after 
cuts, there can be a few of them especially at lower energies) are not taken into account in 
Figs. 10, 11. However, when evaluating $ PT such events are included in the data sample 
2 . All charged particles are used to prepare these plots. No significant differences between 
the distributions of M(p T ) for data and mixed events are observed at any SPS energy, 
indicating the absence of substantial event-by-event fluctuations. The difference between 
the histograms for data and mixed events can be better seen in Fig. 11, which presents 
the ratio of the two. At all energies, the ratio is close to unity, thus confirming that 
average transverse momentum fluctuations are very close to the statistical fluctuations. 
The widths of the distributions of M(p T ) decrease with increasing energy of the colliding 
system as expected from the increasing particle multiplicities. 
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Figure 11: Distributions of mean transverse momentum for data divided by distribution for 
mixed events (see the previous figure). 



The fluctuation measure <3> PT , which is more sensitive to small average pr fluctuations, 
is shown in Fig. 12. The measured values are corrected for limited two-track resolution. 
Results are presented for all charged particles and separately for negatively and positively 
charged particles 3 . For all three charge selections, <3> PT seems to be independent of energy 
and consistent with the hypothesis of independent particle production ($ PT ps 0). The 

2 It was verified that the $ PT value does not depend on whether zero multiplicity events are included or 
not (in agreement with the definition of Q PT ). However, in Table 4, the mean multiplicities of accepted 
particles are shown, where events with zero multiplicity are normally taken into account. Thus, for 
consistency, the code calculating $ PT values includes events with zero multiplicity. 

3 The sample of negatively charged particles is composed mainly of negative pions, whereas the sample 
of positively charged particles is dominated by positive pions and protons. Therefore, the measured 
fluctuations could differ between both charges. Moreover, the sample of all charged particles can include 
additional sources of correlations which arc not present in positively or negatively charged particles, 
separately. Therefore the $ PT measure obtained from all charged particles will not necessarily be a sum 
of $ PT for positively and for negatively charged particles. 
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measured $ PT values do not show any anomalies which might appear when approaching 
the phase boundary or the critical point. However, it should be noted that because of 
the limited acceptance of NA49 and the additional restrictions used for this analysis, 
the sensitivity for fluctuations may be significantly reduced if the underlying scale of the 
correlations is large (see Sec. 2). 
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Figure 12: (Color online) & PT as a function of energy for the 7.2% most central Pb+Pb inter- 
actions. Data points are corrected for limited two-track resolution. Errors are statistical only. 
Systematic errors are given in Table 4. 

Table 4 presents the mean multiplicities of accepted particles, the dispersions <tn = 
y/ (N 2 ) — (N) 2 of the multiplicity distributions, the mean inclusive transverse momenta, 
the dispersions a VT of inclusive transverse momentum distributions, and $ PT values for all 
analyzed data sets. The $ PT values (with their statistical and systematic errors) shown in 
this table have been calculated for all accepted charged particles as well as for negatively 
and positively charged particles, separately. All $ PT values are corrected for the two-track 
resolution effect of the NA49 detector. 
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Energy (GeV) 


(N) 


a N 


(MeV/c) 


a PT (MeV/c) 


$ PT (±stat ±sys) (MeV/c) 


20A (all) 


29 


6 


314 


237 


0.3 ± 0.3 t&g 


20A (-) 


10 


3 


221 


163 


1.3 ± 0.2 


20A (+) 


19 


4 


361 


254 


1.2 ± 0.3 t$ 


30A (all) 


38 


6 


323 


247 


-0.3 ± 0.3 ±H 


30A (-) 


14 


4 


238 


178 


1.1 ± 0.2 


30A (+) 


24 


5 


374 


267 


0.3 ± 0.3 j^g 


40A (all) 


44 


7 


326 


251 


0.6 ± 0.5 ± 0.7 


40A (-) 


17 


4 


242 


182 


1.7 ± 0.4 ± 0.4 


40A (+) 


27 


5 


378 


273 


0.6 ± 0.6 ± 1.1 


80A (all) 


61 


9 


331 


256 


-0.3 ± 0.6 ± 0.9 


80A (-) 


26 


5 


270 


205 


2.4 ± 0.4 ± 0.5 


80A (+) 


35 


6 


378 


280 


-1.4 ± 0.6 ± 0.9 


158A (all) 


77 


10 


331 


255 


0.6 ± 0.7 ± 1.5 


158A (-) 


34 


6 


282 


214 


1.4 ± 0.6 ± 1.0 


158A (+) 


43 


7 


370 


277 


-0.5 ± 0.8 ± 1.7 



Table 4: Measured inclusive and event-by-event parameters for accepted particles (results with 
additional y* cut: see the text for more details). (N), gn, PT, and a PT values are not corrected 
for acceptance. $ PT values are corrected for limited two-track resolution. <E> PT values are given 
with statistical and systematic errors. 



Two-particle correlation plots (for all charged particles) of the cumulant transverse 
momentum variable x are presented in Fig. 13 for 20A, 30A, A0A, 80A, and 158A GeV 
central Pb+Pb collisions. One observes that the plots are not uniformly populated. The 
color scale is the same for all plots in order to check whether the correlation pattern 
changes with the energy. For all five energies, short-range correlations (Bose-Einstein and 
Coulomb effects) are visible as an enhancement of the point density in the region close to 
the diagonal. The analysis of events simulated by the UrQMD model (see below) resulted 
in uniformly populated two-particle correlation plots without diagonal enhancement. 

Figures 14 and 15 present two-particle correlation plots for different charge combi- 
nations. The correlations of both positively and negatively charged particle pairs look 
similar (enhancement along diagonal) and can be explained by the Bose-Einstein effect. 
On the other hand, there is no significant correlation observed for unlike-sign particles. 
The preliminary analysis for (H — ) pairs without azimuthal angle restrictions (see the ac- 
ceptance limits in Fig. 2) showed a small maximum restricted to the region of low x (i.e. 
low transverse momenta) [24]. Such a low x (low px) maximum for unlike-sign particles 
was observed in the STAR data [13] and recently also by the CERES experiment [25]. 
This maximum was interpreted as an effect of Coulomb attraction of particles with differ- 
ent charges and contamination from e + e~ pairs. When the azimuthal angle is restricted 
(Figs. 14 and 15) the maximum at low x is not observed any more. 



17 



J- 1 



80A GeV 



0.8 
0.6 
0.4 
0.2 



I 



1.02 
1.01 
1 

10.99 



* 1 

0.8 
0.6 
0.4 
0.2 



158A GeV 



0.2 0.4 0.6 0.8 1 



i 



1.02 
1.01 
1 

0.99 



0.2 0.4 0.6 0.8 1 



Figure 13: (Color) Two-particle correlation plots (xi,£2) using the cumulant px variable x. 
The bin contents are normalized by dividing with the average number of entries per bin. Plots 
are for all charged particles for central Pb+Pb collisions at 20^4 - 158^4 GeV. Note that the color 
scale is the same in all panels. 
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Figure 14: (Color) Two-particle correlation plots (x\,X2) using the cumulant px variable x. 
The bin contents are normalized by dividing with the average number of entries per bin. Plots 
are for central Pb+Pb collisions at 20^4 GeV for unlike-sign particles (H — ), and for negatively 
and positively charged particles, separately. Note that the color scale can be slightly different 
for different panels. 
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Figure 15: (Color) Two-particle correlation plots (xi,X2) using the cumulant px variable x. 
The bin contents are normalized by dividing with the average number of entries per bin. Plots 
are for central Pb+Pb collisions at 158A GeV for unlike-sign particles (H — ), and for negatively 
and positively charged particles, separately. Note that the color scale can be slightly different 
for different panels. 
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5.2 Comparison with the UrQMD model 

The measured $ PT values have been compared with the predictions of the ultrarelativis- 
tic quantum molecular dynamics (UrQMD) model [26, 27], a transport model producing 
hadrons via formation, decay, and rescattering of resonances and strings. In the UrQMD 
model, no fluctuations due to a phase transition are incorporated. However, resonance 
decays and effects of correlated particle production due to energy (momentum) and quan- 
tum number conservation laws are included. In the analysis default parameters of the 
UrQMD model were used (meson- meson and meson-baryon scattering included). For 
each energy the most central 7.2% interactions were selected, in accordance with the real 
NA49 events. 

In the study of the NA49 data, the $ PT measure was calculated from all charged 
particles, consistent with originating from the main vertex. This means that mostly main 
vertex pions, protons, kaons, and their antiparticles are used in the analysis, because 
particles coming from the decays of Kg, A, H and Q are suppressed by the track selection 
cuts. Therefore, also in the analysis of the UrQMD events only charged pions, protons and 
kaons were taken into account for evaluating <3> PT . Due to the specific parameter applied 
in the UrQMD model - the time during which all particles are tracked (taken as 80 fm/c) 
- the list of generated kaons, pions and (anti)protons does not contain the products of 
weak decays. In the analysis of the UrQMD events, the same kinematic restrictions were 
applied as in the case of the NA49 data. In particular, the selected (4>,Pt) acceptance is 
common for all five energies, and as a result of all the kinematic and acceptance cuts, only 
about 5-6% of all charged particles are used to study transverse momentum fluctuations. 
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Figure 16: (Color online) Comparison of & PT as a function of energy from data (data points, 
corrected for limited two-track resolution) with UrQMD model calculations (black lines) with 
acceptance restrictions as for the data. The panels represent results for all charged (left), 
negatively charged (center) and positively charged particles (right). 



Figure 16 compares the the energy dependence of $ PT for the data and for the UrQMD 
model. Similar to the data, the UrQMD model does not show any significant energy de- 
pendence of transverse momentum fluctuations for any charge selection. The range of 
$ PT values obtained from the UrQMD calculations is similar to that found in the data. 
Quantitative comparisons are not conclusive, since in spite of the same kinematic and 
acceptance restrictions the model does not contain the effects of Bose-Einstein correla- 
tions and Coulomb interactions. Correlations implemented in the UrQMD model (due 
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to conservation laws and hadron resonance production and decays) apparently do not 
increase the observed $ PT values. In particular the observed differences between model 
and experimental data do not point to any onset of unusual fluctuations as expected for 
hadronic freeze-out close to the critical point. 



The sensitivity of the results on the energy dependence of $ PT to variations of the upper 
Pt cut was checked by repeating the analysis for three different upper p T cuts. Figures 17, 
18 and 19 show the dependence of <3> PT on energy for upper p T cuts of 750 MeV/c, 500 
MeV/c and 250 MeV/c, respectively. The measurements (data points) are corrected for 
limited two-track resolution. The lines represent the UrQMD predictions with the same 
kinematic and acceptance restrictions. The results with a decreased upper p? cut do not 
show any significant energy dependence, which is similar to the corresponding results for 
the wide-pr interval (pt < 1500 MeV/c). The UrQMD results seem to lie systematically 
below the NA49 data. This effect might be explained by the neglect of Bose-Einstein and 
Coulomb correlations (see also Fig. 13) in the model. 
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Figure 17: (Color online) <& PT as a function of energy for the 7.2% most central Pb+Pb inter- 
actions, with additional cut px < 750 MeV/c. Data (points) are corrected for limited two-track 
resolution. Lines show results of UrQMD calculations with the same kinematic and acceptance 
restrictions. Errors are statistical only. 



5.3 Search for the critical point 

Recent results on the energy dependence of hadron production properties at the CERN 
SPS [2, 3] indicate that deconfined matter is formed at the early stage of central Pb+Pb 
collisions at energies as low as 30A GeV. Thus, at higher collision energies, the expanding 
matter crosses the phase boundary between deconfined matter and hadron gas and may 
freeze-out close to it. The nature of the transition is expected to change with increasing 
baryochemical potential At high potential, the transition is believed to be of the 
first order with the end point of the first order transition line being a critical point 
of the second order. A characteristic property of the second-order phase transition is 
a divergence of the susceptibilities. Consequently, an important signal of the critical 
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Figure 18: (Color online) Same as Fig. 17, but with additional cut pt < 500 MeV/c. 
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Figure 19: (Color online) Same as Fig. 17, but with additional cut pt < 250 MeV/c. 



point are large fluctuations, in particular, an enhancement of fluctuations of transverse 
momentum and multiplicity [5]. These predictions are to a large extent qualitative, as 
QCD at finite temperature and baryon number is one of the least accessible domains of 
the theory. Thus quantitative estimates of the critical point signals were performed within 
QCD-inspired models [5, 28, 29]. 

It was found that assuming a correlation length £ = 6 fm and freeze-out at the critical 
point $ PT in full phase space may be increased by about 40 and 20 MeV/c for all charged 
and for like-charge hadrons, respectively [5] . These predictions are independent of particle 
charge and scale with multiplicity. Therefore $ PT for all charged hadrons is two times 
larger than for like-charge hadrons. The value of $ PT is expected to decrease because of the 
limited acceptance, namely, by a factor of 0.6 and 0.4 for the NA49 acceptance in rapidity 
and in azimuthal angle, respectively. The particle correlator [5] at the critical point was 
calculated to be independent of azimuthal angle and charge of hadrons, whereas it changes 
with rapidity. Therefore the predicted $ PT values scale with azimuthal acceptance 
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in this analysis, independent of energy (see Fig. 2)). Together, these reduce the expected 
signal to about 10 and 5 MeV/c for all charged and like-charge hadrons, respectively 
[5] (private communication). These estimates are based on a correlation length £ = 6 
fm. However, because of the finite lifetime of the fireball, the correlation length may not 
exceed 3 fm/c [30]. Since the fluctuations scale with the square of the correlation length 
the expected signal is further reduced by (6/3) 2 = 4. This leads to estimates for <3> PT of 
10 MeV/c for all charged hadrons within full phase space and about 2.4 MeV/c for all 
charged hadrons and 1.2 MeV/c for like-sign charged hadrons in the acceptance of the 
analysis. Theoretical estimates based on lattice QCD calculations locate the critical point 
at T « 162 MeV and \i B « 360 MeV [4]. Guided by the considerations of Ref. [31], we 
parametrize the increase of $ PT due to the critical point by Gaussian shapes in T and \xb 
with a{T) « 10 MeV and a(fi B ) ~ 30 MeV, respectively. 
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Figure 20: (Color online) <3? PT as a function of baryochemical potential (jj,b values from statis- 
tical hadron gas model fits [32]) for the 7.2% most central Pb+Pb interactions. Results with 
Pt < 1.5 GeV/c (as in Fig. 12). Data (points) are corrected for limited two-track resolution. 
Errors are statistical only. The Gaussian curves show estimated $ PT values in the case of the 
existence of the critical point (see text for details). 



The results on $ PT in central Pb+Pb collisions at the SPS energies for 0.005 < pt < 
1.5GeV/c and the restricted acceptance of the analysis are shown as a function of \x& in 
Fig. 20. The baryochemical potential ^ at each energy was obtained from fits of the 
statistical hadron gas model to the measured particle yields of NA49 [32]. The predictions 
for the effects of the critical point for the two assumptions on the correlation length (3 
and 6 fm) are illustrated by the dashed and full Gaussian curves, respectively. The base 
lines of the model curves are set to the mean level of $ PT in the SPS energy range (here 
the mean level is calculated separately for all charged and negatively charged particles). 
The data do not support the predictions for the critical point. Neither a nonmonotonic 
increase of $ PT nor a characteristic difference between $ PT for all charged and like-charge 
hadrons are observed. However, the freeze-out in central Pb+Pb collisions appears [32] to 
take place at temperatures significantly below the transition temperature and thus could 
make the critical point signal invisible. This possibility motivates the continuation of 
the search for the critical point in NA61/SHINE [33] by a scan in both temperature and 
baryochemical potential which will be performed by changing the collision energy and size 
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of the colliding nuclei. 

Since publication of the NA49 results on multiplicity fluctuations [6], the estimates for 
the effects of the critical point have been updated. In the following, we therefore present 
these new predictions for multiplicity fluctuations for consistency. The estimates of the 
critical point signals were again performed [5] (private communication) assuming first a 
correlation length £ = 6 fm, and they suggest an increase of the scaled variance of the 
multiplicity distribution u in full phase space by 2 for all charged and by 1 for same-charge 
particles (uncorrelated particle production results in u = 1). The limited acceptance used 
in the multiplicity fluctuation analysis [6] leads to correction factors for the increase of 
0.6 and 0.7 due to the rapidity and azimuthal acceptance, respectively. In the analysis of 
multiplicity fluctuations, the azimuthal acceptance was chosen dependent on the energy 
with a mean of about 70%. Thus the measured multiplicity fluctuations are expected 
to be increased by 0.84 and 0.42 for all charged and like-charge hadrons, respectively. 
For a correlation length £ = 3 fm, the increase should be four times smaller, i.e., 0.5 for 
all charged and 0.25 for same-charge particles in full phase space and, respectively, 0.21 
and 0.10 in the NA49 acceptance. Figure 21 shows the experimental results for u as a 
function of baryochemical potential in the rapidity interval 1 < y* < yl eam for the 1% 
most central Pb+Pb collisions [6]. The expectations for the critical point are represented 
by the Gaussian curves above base lines at the level of the mean measured values of u. 
Similar to transverse momentum fluctuations, multiplicity fluctuations do not exhibit the 
effects expected for hadronic freeze-out close to the critical point. However, one should 
note that the predicted width of the critical point signal is comparable to the size of the 
fiB steps in the NA49 energy scan. In future experiments, one should consider scanning 
with narrower steps in \lb- 



1.5 



fluct. at CEP: 


• all charged 




fluct. at CEP: 


■ negatively charged 


all charged; %=Z fm 

— all charged; ^=6 fm / 

* 
i 

— i 


1 < y < y 

-» ti beam 
* 

t 

\ 




neg. charged; ^=3 fm 
neg. charged; ^=6 fm 


▲ positively charged 
1 < y < y 

% beam 




..^^T^r. __ __ _. 




* 

* 






i — • — i i-m-i 








' 1 1 


1 1 


' 1 


1 


i i 


200 300 


400 500 200 


300 


400 500 




M-g [MeV] 






lx B [MeV] 



Figure 21: (Color online) uj [6] as a function of baryochemical potential (fiB values from statis- 
tical hadron gas model fits [32]) for the 1% most central Pb+Pb interactions in 1 < y* < y^ eam - 
Errors are statistical only. The Gaussian curves show estimated u values in the case of the 
existence of the critical point (see text for details). 
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5.4 Comparison with other experiments 



Event-by-event transverse momentum fluctuations have been studied by other experi- 
ments both at SPS and at RHIC energies. This section compares the NA49 results with 
those from the CERES and the STAR experiments. In such comparisons, one has to 
remember that the magnitude of any fluctuation measure may depend on the acceptance 
used in the analysis. The CERES experiment measured somewhat higher values of pr 
fluctuations in the midrapidity region for central Pb+Au collisions at 40A, 80v4, and 158A 
GeV. However, no significant energy dependence was found for either the E PT (%) or the 
$ PT fluctuation measures [16, 34]. Figure 22 compares NA49 and CERES [16] results on 
the energy dependence of the $ PT measure. One observes only very weak (if any) energy 
dependence of $ PT over the whole SPS energy range. It should, however, be stressed that 
quantitative comparison of $ PT values from NA49 and CERES is hampered by different 
acceptances (NA49: forward rapidity and limited azimuthal angle, CERES: midrapidity 
and complete azimuthal acceptance). 
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Figure 22: (Color online) <& PT as a function of energy measured for all charged particles by the 
NA49 and CERES experiments. The NA49 points are obtained for the forward-rapidity region 
in a limited azimuthal angle acceptance; the CERES data [16] are calculated for the midrapidity 
region within a complete azimuthal acceptance. All data points include small-scale correlations 
(HBT, etc.). 

The <3> PT measure is close to zero at SPS energies, but the STAR [35] results show 
a strong increase of $ PT from top SPS to RHIC energies. The energy dependence of 
$ PT is shown for central Pb+Pb (Pb+Au, Au+Au) collisions in Fig. 23. The data of 
NA49 (this paper), CERES, [16] and STAR [35] are compiled. The STAR results are 
presented for the related measure Aa PT . n instead of $ PT . The difference between both 
measures is smaller than a few percent for high multiplicity collisions as considered here 
(no significant difference within the SPS energy domain and only few-percent differences 
at RHIC energies [35]). 

The results of NA49 and CERES at v /saT/v=17.3 GeV differ significantly from those of 
STAR at v /J/v r /v=19.6 GeV. These differences are likely due to differences in the acceptance 
in momentum space. This is illustrated in Fig. 24, where $ PT (A<r PT:n ) is plotted as a 
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Fi gure 23; (Color online) ^p T or Acjp T:n as a function of energy measured for all charged 
particles by the NA49, CERES, and STAR experiments. The NA49 points are obtained for the 
forward-rapidity region in a limited azimuthal angle acceptance; the CERES [16] and STAR (Fig. 
3, left, in Ref. [35]) points are calculated for midrapidity regions within complete azimuthal 
acceptances. All data points include small-scale correlations (HBT, etc.). STAR data are 
corrected for tracking inefficiency 



function of the fraction of accepted charged particles 4 . The values of <& PT (Aa PT . n ) 
increase with increasing fraction of accepted particles. 

The striking feature of the data presented in Fig. 23 is the rapid increase of $ PT 
(A(Tp T:n ) with increasing collision energy. For a given experiment, the fraction of ac- 
cepted charged particles is constant (NA49), increases (CERES), or decreases (STAR) 
with increasing collision energy 5 . It was already explained in Sec. 2 that $ PT is indepen- 
dent of acceptance changes provided the correlation scale is much smaller than the size of 
the acceptance region. On the other hand, when the range of correlations is significantly 
larger than the acceptance, $ PT linearly decreases with decreasing fraction of accepted 
particles. The increase of $ PT (Acr PT:n ) with energy is not caused by an increasing mag- 
nitude of short range (compared to the acceptance) correlations. As mentioned above, 
these correlations, if dominant, would lead to an independence of <3> PT (A<7 PT:n ) of the 
acceptance, in contrast to the results shown in Fig. 24. Curiously enough, the increase of 

4 The plot shows A+A data at the top SPS energy and the lowest RHIC energy. The total multiplicities 
of charged particles (N c h) have been obtained from the UrQMD model with centralities corresponding 
to the NA49, the CERES, and the STAR data. The following numbers were estimated: N ch = 1665 for 
5% most central Pb+Pb (NA49), N ch = 1567 for 7.2% most central Pb+Pb (NA49), N ch = 1594 for 5% 
most central Pb+Au (CERES), N ch = 1575 for 6.5% most central Pb+Au (CERES), and N ch = 1401 for 
20% most central Au+Au (STAR). The values of $ PT (Aer pr:n ) and multiplicities of accepted particles 
were taken from [10, 16, 35] and from this paper. 

5 The fraction of accepted charged particles in NA49 is about 5%, in CERES it increases from about 
5.8% at 40A GeV to about 10% at 80A and 158A GeV, and for STAR this fraction decreases from about 
31% at the lowest RHIC energy to approximately 17-18% at the top RHIC energy. In order to obtain the 
above numbers the total multiplicities of charged particles were generated within the UrQMD model with 
centralities corresponding to that of the NA49, the CERES, and the STAR data presented in Fig. 23. 
The multiplicities of accepted charged particles were taken from [16], [35] and from this analysis. 
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Figure 24: (Color online) <J? PT or Aa PT - n as a function of percent of charged particles accepted in 
the analysis. Data are presented for the top SPS and the lowest RHIC energies. The NA49 points 
(Ref. [10] and this analysis) are obtained for the forward-rapidity region in a limited azimuthal 
angle acceptance; the CERES [16] and STAR [35] points are calculated for midrapidity regions 
within complete azimuthal acceptances. All data points include small-scale correlations (HBT, 
etc.). 



Acr pT:ri with growing ^snn for the STAR data is associated with a decrease of the fraction 
of accepted particles. Consequently, we conclude that the increase is caused by a growing 
magnitude of medium- and long-range correlations (correlation length comparable to or 
larger than the acceptance). This effect may result from an increasing contribution of 
particles originating from decays of high mass resonances and/or from (mini)jet fragmen- 
tation [35]. 

In contrast to the $ PT (Aa PT;n ) fluctuation measure, £ PT used by the CERES exper- 
iment, does not show any dramatic differences when going from SPS to RHIC energies 
[16, 25] (see also Fig. 25 6 ). Also the PHENIX results confirm that the magnitude of the 
E PT (%) fluctuation measure exhibits only a small variation in a very wide energy range 
between -^snn = 22 and 200 GeV [37]. It appears that the increase of $ PT (A<7 PT:n ) with 
the accepted particle multiplicity is approximately linear (Fig. 24). Consequently, S PT (%), 



which is proportional to y$ PT /(A^), becomes approximately independent of multiplicity 
or collision energy. 

A direct comparison of event-by-event px fluctuations between different experiments, 
based on the method of determining the magnitude of residual temperature fluctuations 



6 The NA49 £ Pr points were estimated from the values of <J> Pr using the following equations: 

a l T ,dyn - ^"(n)^ and finall y s pt = s 9n{pl T dyn ) ■ ^ dyn [16]. For a given A$ PT value the errors 
of £ pr (AT, PT upper and A£ pr /otuer ) were calculated as: AZ PT upper = |S Pr ($ PT + A$ Pr ) - S Pr ($ PT )| 
and AY, PT lower = |S PT ($ PT - A$ PT ) - £ Pt ($ Pt )|. Finally the NA49 $ PT ± A$ PT corresponds to 

+ A£ upper 

£ PT _ AS PT ioroer . The same procedure was applied for statistical and systematic errors. Although the 
statistical and systematic errors of $ PT are comparable at CERES and NA49, the errors of E PT are much 
higher in NA49. This is due to the fact that S PT is proportional to sgn(Q PT ) ■ yJ\$ PT | and for <I> PT close 
to zero (the case of NA49) a small variation of $ PT results in much higher changes of E PT . 
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Figure 25: (Color online) E PT (%) as a function of energy measured for all charged particles 
by the NA49, CERES, and STAR experiments. The NA49 points are obtained for the forward- 
rapidity region in a limited azimuthal angle acceptance; the CERES ([16]; point for 8% most 
central Pb+Au collisions taken from Ref. [25]) and STAR [36] points are calculated for midra- 
pidity regions within complete azimuthal acceptances. 

[38], has been presented in Ref. [39]. For the most central interactions, the estimated 
temperature fluctuations are on the level of 1.8%, 1.7%, 1.3%, and 0.6% for PHENIX 
(y^^ = 200 GeV), STAR (y/s]^ = 130 GeV), CERES (^/W = 17 GeV), and NA49 
(y/s NN = 17 GeV), respectively. This analysis confirms that, indeed, there is no significant 
energy dependence of temperature fluctuations between SPS and RH1C energies. 



6 Summary and outlook 

Transverse momentum event-by-event fluctuations were studied for central Pb+Pb in- 
teractions at 20A, 30A, 40A, 80A, and 158A GeV. The analysis was limited to the 
forward-rapidity region (1.1 < y* < 2.6). Three fluctuation observables were studied: 
the fluctuations of average transverse momentum (M(pp)) of the event, the $ PT fluctua- 
tion measure, and transverse momentum two-particle correlations. The following results 
were obtained: 

• The distribution of event-by-event average transverse momentum M(pj) is close to 
that for mixed events, indicating that fluctuations are predominantly of statistical 
nature. 

• The fluctuation measure $ PT has values close to zero and shows no significant energy 
dependence. 

• Two-particle transverse momentum correlations are found to be negligible except 
for contributions from Bose-Einstein correlations. 

• The UrQMD model reproduces the trend of the data. 
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• The observed smallness of pt fluctuations in the SPS energy range is consistent with 
other published measurements. 

• No sudden increase or nonmonotonic behavior was observed in the energy depen- 
dence of $ PT nor the scaled variance uo of multiplicity fluctuations. Such effects have 
been suggested for freeze-out near the critical point of QCD. 

Since the limited acceptance of NA49 is expected to reduce the signal from possi- 
ble critical fluctuations, measurements in larger acceptance are still needed. Moreover, 
employing intermediate size nuclei may possibly move the freeze-out point closer to the 
critical point. These considerations motivated several new experiments planned in the 
SPS energy region: NA61 at the CERN SPS [40, 41, 33, 42], STAR and PHEN1X at 
the BNL RHIC [43], MPD at the Joint Institute for Nuclear Research (JINR) Nuclotron- 
based Ion Collider Facility (NICA) [44] , and CBM at the GSI Facility for Antiproton and 
Ion Research (FAIR) heavy-ion synchrotron (SIS-300) [45]. The future results, together 
with the existing data, will cover a broad range in the (T, plane and should lead to 
significant progress in the search for the critical point. 
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